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Synthetic biology holds promise as both a framework for rationally engineering biological systems 
and a way to revolutionize how we fundamentally understand them. Essential to realizing this 
promise is the development of strategies and tools to reliably and predictably control and charac- 
terize sophisticated patterns of gene expression. Here we review the role that RNA can play towards 
this goal and make a case for why this versatile, designable, and increasingly characterizable mole- 
cule is one of the most powerful substrates for engineering gene expression at our disposal. We dis- 
cuss current natural and synthetic RNA regulators of gene expression acting at key points of con- 
trol - transcription, mRNA degradation, and translation. We also consider RNA structural probing 
and computational RNA structure predication tools as a way to study RNA structure and ultimate- 
ly function. Finally, we discuss how next-generation sequencing methods are being applied to the 
study of RNA and to the characterization of RNA's many properties throughout the cell. 
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1 Introduction 

The field of synthetic biology is emerging to solve some of 
the fundamental challenges that arise when we try to har- 
ness the powerful, yet extremely complex, capabilities of 
natural biomolecular systems. The applications that syn- 
thetic biologists are pursuing are diverse and range from 
optimizing the biosynthetic capabilities of cells to pro- 
duce fuels and drugs, to harnessing the ability of cells to 
act as environmental biosensors, to manipulating cellular 
communities and tissues for medical purposes [1, 2]. 
While the challenges raised by these applications 
demand a sophistication of engineering on many levels, 
at the heart of all of these applications is the precise con- 
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trol of gene expression - turning the right genes on at the 
right times and at the right levels. Because of this central 
role of gene expression, one of the grand challenges of 
synthetic biology has become the development of an 
engineering framework that allows the design, construc- 
tion, and characterization of genetic networks that reli- 
ably implement specified genetic control. 

While great strides have been made towards this goal 
by focusing on protein systems (reviewed in [2, 3]), there 
is a case to be made for the power of RNA-based gene reg- 
ulators as fundamental components of the synthetic biol- 
ogy toolbox. To start with, RNA's role as a versatile regu- 
lator of cellular genomes is being uncovered at an accel- 
erating pace [4], and there are now well studied RNA 
mechanisms that regulate almost all aspects of gene 
expression [5-9]. We also have a rapidly emerging bio- 
physical understanding of the relationship between RNA 
structure and function that has served as an excellent 
starting point for computational RNA design [10-14]. But 
perhaps even more compelling is our ability to character- 
ize RNAs in unprecedented breadth and scale with the 
advent of next-generation sequencing (NGS). In fact, 
there has been an explosion of sequencing-based tools 
that can characterize almost all aspects of RNAs - from 
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cellular abundances, to locations of protein and ribosome 
binding, to structures and interactions [15-24] - all with 
variations on an experiment that is rapidly becoming 
cheaper and faster [25]. 

We therefore hypothesize that RNA represents one of 
the most powerful substrates at our disposal with which 
to develop a discipline of engineering gene expression 
inside cells, because it is versatile, designable, and easily 
characterizable. It is our goal in this review to present the 
case for RNA by focusing on these three features in turn. 
We start with a survey of what is known about natural 
RNA regulatory mechanisms and an update on their lat- 
est engineering. We then give a brief introduction and 
overview of computational and experimental RNA struc- 
ture determination techniques and how they fit into the 
RNA engineering design paradigm. Finally, we provide an 
introduction to several important NGS characterization 
techniques for RNA, with pointers to resources for start- 
ing these experiments, and thoughts on how they can be 
applied to RNA engineering problems. To make our dis- 
cussion concrete, we focus on prokaryotic RNA-based 
regulation, although much progress has been made in 
engineering RNA regulation in eukaryotes [26-32]. 

2 Natural and engineered RNA mechanisms 
for control of gene expression 

Sixty years of molecular biology has taught us that gene 
expression can be thought of as a series of biochemical 
conversions from DNA to RNA and in many cases to pro- 
tein (Fig. 1). While this is a simplified model, it is a useful 



framework for establishing the control points of basic 
gene expression. In prokaryotes, the most obvious points 
of control are RNA synthesis (transcription), RNA degra- 
dation, protein synthesis (translation), and protein degra- 
dation. Choosing different points of control can have dif- 
ferent ramifications for the dynamics and steady-state 
levels of gene expression [33]. There are also a variety of 
different mechanisms that can regulate each point of con- 
trol. These factors, combined with potential issues of 
compatibility between different mechanisms, create a 
rich optimization challenge for the biomolecular engineer 
wishing to precisely control gene expression. While pro- 
teins play fundamental roles in gene expression, RNAs 
have also been found to regulate multiple points of con- 
trol, namely transcription, RNA degradation, and transla- 
tion. For each control point, we discuss the main natural 
classes and mechanisms of RNA-based gene regulation 
and those that have been recently engineered by syn- 
thetic biologists (Fig. 2A and B). 

2.1 Translational control of gene expression by RNA 

2. 7. 7 General roles of RNA structure in translation 

Translation initiation has long been known to be a rate- 
limiting step of gene expression [34, 35]. Translation is ini- 
tiated through an interaction between the 16S ribosomal 
RNA (rRNA) and a ribosome binding site (RBS) that lies in 
the 5' untranslated region (5' UTR) of a prokaryotic mes- 
senger RNA (mRNA). Evidence suggests that the ribo- 
some makes contacts within a -52 nucleotide mRNA 
region [36], within which is an eight nucleotide Shine- 
Dalgarno sequence (SD) that is partially complementary 



Trans: 
sRNA 




1 



Cis: 

Riboswitch 
Ribozyme 
Leader peptide 



Promoter 



RBS 



mRNA 
Degradation 



Trans: 
sRNA 




11 





Protein 
Degradation 



0 



Trans: 
sRNA 
Ribozyme 



0 



Figure 1. The central synthesis and degradation steps of prokaryotic gene expression. The points of control where cis- and trans-acting RNA regulators 
exert their regulation are highlighted. Pointed arrows indicate activation while blunt arrows indicate repression. Dashed arrow indicates a gap in regulatory 
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to the 3' end of the 16S rRNA [37, 38]. It is well established 
that variations of the SD sequence away from a comple- 
mentary match to the 16S 3' end have an effect on trans- 
lation rate. Moreover, it has also long been known that 
RNA secondary structures, defined by the pattern of 
canonical base pairs (A pairs to U, G pairs to C or U), in the 
region surrounding the RBS and the start codon can alter 
RBS accessibility and thus tune translation initiation rate 
[39]. Recent evidence by Kudla et al. [40] suggests that 
this may be a general principle of translation regulation. 
In this work, the authors showed that for a library of 154 
synonymously coded green fluorescent protein (GFP) 
variants, the strongest positive correlate of high GFP 
expression was weaker predicted RNA structure near the 
RBS and start codon [40]. 

More recently, synthetic biologists have sought strate- 
gies to minimize the effect of 5' UTR structures. Mutalik 
et al. [41] recently demonstrated one such strategy that 
uses two SD sequences in the mRNA - one that controls 
translation of a leader peptide, and a second that controls 
translation of the downstream gene. Co-translation of the 
leader peptide disrupts secondary structures within the 
downstream RBS, thus minimizing secondary structure 
effects that arise from different combinations of genetic 
elements. Alternative strategies have sought to physical- 
ly remove excess 5' UTR sequences of mRNAs by cleav- 
age via ribozymes [42] , or with ribonucleases derived from 
the bacterial clustered regulatory interspaced short palin- 
dromic repeats (CRISPR) defense system [43]. Moreover, 
as discussed in Section 3, computational tools are now 
being developed to design 5' UTR sequences with RNA 
folding in mind so that translation initiation rates can be 
predictably tuned [10]. 

While the presence of general RNA structure in the 5' 
UTR can adjust RBS accessibility, many natural and syn- 
thetic RNA regulators of translation specifically target the 
RBS and start codon with either trails-acting RNA-RNA 
interactions or cis-acting RNA structures that can change 
conformations and thus switch translation on or off, as 
described below. 

2. 7.2 Translational control: Small RNAs 

Small RNAs (sRNAs) are a diverse and abundant class of 
short (50-100 nucleotides) transcripts that directly target 
specific mRNAs through intermolecular base-pairing 
interactions [6, 44] (Fig. 2A). They can either be cis- 
encoded antisense RNAs, and thus fully complementary 
to their targets, or trans-encoded at a distant location 
with only partial target complementarity [5]. In the latter 
case, these sRNAs often require the RNA chaperone pro- 
tein Hfq hypothesized to stimulate base-pairing by 
recruiting both sRNA and target mRNA to its surface, 
induce structural rearrangements in either the sRNA or 
mRNA, stabilize bound sRNA from ribonuclease degra- 
dation and stimulate degradation of sRNA and mRNA 
hybrids via recruitment of ribonucleases [45] . 



By far the most abundant natural mechanism of gene 
regulation by sRNA is at the post-transcriptional level. 
The simplest mechanism involves sRNA binding and 
occlusion of the RBS or start codon of a target mRNA [46] 
(Fig. 2A), however binding outside of this region can also 
alter translation [47-48]. Figure 2A highlights several 
sRNA translational regulators and the natural genes they 
target such as IS10 RNA-OUT/RNA-IN [49], and sok/hok 
[50]. 

In addition to acting as repressors, in a few cases 
sRNAs have been shown to be translational activators 
[51]. Generally known as the anti-antisense mechanism, 
translational activation occurs when sRNA-mRNA base- 
pairing interactions induce the formation of alternative 
mRNA structures that expose the RBS. Notable examples 
include the sRNA RNAIII activation of hla in Staphylo- 
coccus aureus [52] and DsiA activation of ipoS in E. coli 
[53] (Fig. 2A). 

One of the earliest strategies of engineering synthetic 
sRNA-like regulation was the development of riboregula- 
tors [54] (Fig. 2A). Riboregulators are in essence a cleaner 
implementation of the natural sRNA regulation principle. 
Target mRNAs were designed to be cis-repressed RNA 
(crRNA) by the inclusion of a 5' intrinsic stem-loop struc- 
ture that encompassed the RBS. Activation of translation 
was achieved by the introduction of an sRNA analog 
called a trans-activating RNA (taRNA) designed to target 
and hybridize this stem-loop to expose the RBS. More- 
over, these riboregulators have been further engineered to 
provide orthogonal crRNA and taRNA pairs that have 
been effectively combined in higher-order devices such 
as a four input genetic switchboard [55] and a synthetic 
counter [56]. This approach was recently extended by 
Rodrigo et al. [12], who developed an algorithm to com- 
putationally search for additional crRNA/taRNA pairs 
using a structure-guided design. 

Additional recent strategies have been to create vari- 
ants of natural sRNAs. For example, Mutalik et al. [57] 
engineered the copy number control element from the 
transposon IS10, consisting of the sRNA (RNA-OUT) and 
the target mRNA (RNA-IN; Fig. 2A). In a mechanistically 
driven approach, RNA-IN and RNA-OUT variants were 
generated by mutating the sequences known to govern 
the specificity of the interaction. By gathering gene 
expression data for a 52 x 52 matrix of different RNA- 
IN/RNA-OUT pairs, the authors were able to construct a 
statistical model to find sequence elements and predict- 
ed structural features that were critical to function. In par- 
ticular, they found that the overall hybridization energy of 
the RNA-IN/RNA-OUT pair, and a 5 bp region that seeds 
the hybridization, were key features that correlated with 
function [57]. 

Synthetic Hfq-binding sRNAs have also been designed 
by using existing sRNAs as scaffolds to provide Hfq bind- 
ing sites and transcriptional terminators while modifying 
the mRNA target binding sequence. Sharma et al. [58] 
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screened libraries of randomized target binding sequence 
to generate sRNAs able to repress a variety of endogenous 
E. coli genes. Alternatively, Na et al. [59] rationally 
designed target sequences complementary to either the 
translation initiation or coding regions of target mRNAs. 
These synthetic sRNAs were used in a metabolic engi- 
neering approach to combinatorially repress the expres- 
sion of chromosomal genes and identify E. coli strains 
with the highest production of a desired metabolite. One 
of the major advantages of using sRNA for the repression 
of endogenous genes is that no chromosomal modifica- 
tion is required to generate knock-down strains [59]. 

Going further, natural sRNAs have also recently been 
engineered to change structure and function in response 
to external ligands (Fig. 2A). Qi et al. [60] fused the IS10 
RNA-OUT hairpin to a well characterized RNA aptamer - 
an RNA structure that can bind to a specific molecule. In 
this design, the loop of the aptamer hairpin was made 
complementary to the RNA-OUT loop to form a pseudo- 
knot interaction in the absence of ligand, rendering the 
antisense non-functional. Ligand binding to the aptamer 
prevented pseudoknot formation, thus allowing the RNA- 
OUT to fold properly and repress its target. This trans-act- 
ing, ligand-dependent conditional folding strategy was 



actually inspired by natural ra's-acting RNA regulators 
called riboswitches, which we discuss next. 

2.1.3 Translational control: Riboswitches 

Riboswitches are a class of RNA regulators that control 
gene expression in response to ligands, and appear to be 
used ubiquitously across bacteria (Fig. 2A) [61]. Simple 
riboswitches are composed of two domains: an aptamer 
and a gene expression platform [62]. These domains are 
structurally linked, with ligand-dependent rearrange- 
ments in the aptamer domain causing RNA conforma- 
tional changes within the expression platform that lead to 
different gene expression outcomes [62]. Riboswitches 
that regulate translation have expression platforms that 
fold to either occlude or expose the RBS, and can be either 
activators such as the S-adenosylhomocysteine (SAH) 
riboswitch ahcY [63] (Fig. 2A), or repressors as in the case 
of the adenosylcobalam [64] and thiamine [65] ribo- 
switches (Fig. 2A). Some riboswitches such as the glycine 
riboswitch have tandem copies of the same aptamer 
before the expression platform, allowing cooperative 
response to ligands [66]. There are also examples of 
riboswitches with two different aptamers that allow inte- 
gration of multiple inputs as in the metE riboswitch [67]. 
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As suggested by their ubiquity and utility in nature, 
synthetic riboswitches hold great potential for engineer- 
ing biological systems that can sense and respond to envi- 
ronmental or intracellular inputs. Interestingly, natural 
riboswitches were only discovered after they were first 
engineered in the laboratory by combining an aptamer 
that binds to ATP and a self-cleaving ribozyme [68]. Since 
then, progress has been made in developing systematic 
selection and engineering strategies to fuse different 
aptamers to appropriate expression platforms to make 
synthetic riboswitches (Fig. 2A). For example, Desai and 
Gallivan [69] screened a library consisting of the theo- 
phylline aptamer upstream of a variable linker sequence 
and an RBS to find a theophylline riboswitch with 36-fold 
activation. In another example, Suess et al. [70] rationally 
designed a theophylline riboswitch that had alternative 
stem-loop folds in close proximity to the RBS. Further 
progress has been aided by the development of high- 
throughput functional screening assays to rapidly screen 
sequence variants including colorimetric [71] and flow- 
cytometry [72] based screens, and even selection strate- 
gies based on cell motility [73]. These approaches should 
aid in the development of riboswitches from aptamers dis- 
covered from the systematic evolution of ligands by expo- 
nential enrichment (SELEX) method [74-76], an in vitro 
technique that has allowed the discovery of aptamers that 
bind to a wide array of ligands [76]. In fact, riboswitches 
have already been constructed from a number of different 
aptamers including those responsive to antibiotics [77, 
78], proteins [79] and other small molecules [68, 78, 80]. 

2.7.4 Translational control: Ribozymes 

Ribozymes are naturally occurring RNA sequences that 
have the ability to catalyze the cleavage or ligation of RNA 
strands (for review see [81]). Most synthetic ribozymes 
have utilized the hammer-head ribozyme, which has also 
been fused with aptamers to create ligand-controlled 
ribozymes called aptazymes [68]. For example, Wieland 
and Hartig [82] developed a theophylline-responsive 
aptazyme that cleaved structures occluding the RBS, 
thereby activating gene expression (Fig. 2A). More 
recently, ribozymes have been developed that self-cleave 
in response to trans-acting sRNAs that are functional in 
vivo, providing the potential for coupling ribozymes to 
existing and synthetic sRNAs for the development of RNA 
networks [83] (Fig. 2A). 

2.2 Control of gene expression 
by mRNA degradation 

2.2. 7 General roles of RNA structure 
in mRNA degradation 

In bacteria, mRNA degradation is primarily achieved by 
the RNA degradosome, a protein complex composed of 
RNase E, polynucleotide phosphorylase (PNPase), and 
RhlB RNA helicase [84]. Degradation is a multiple step 



process, with initial RNase E cleavage being the most rate 
limiting. Interestingly, since RNase E cleaves single- 
stranded RNA, secondary structures can inhibit its activ- 
ity. Indeed, this has been effectively used to engineer spe- 
cific secondary structures within the 5' UTR of mRNAs to 
increase stability, and therefore overall gene expression, 
over a 10-fold range [85]. 

2.2.2 mRNA degradation control: sRNA 

Natural sRNAs that regulate translation also affect mRNA 
degradation. In fact, sRNA-mRNA interaction is typically 
followed by rapid degradation by RNase E or RNase III [5]. 
Moreover, Hfq has been shown to interact with compo- 
nents of the degradosome, highlighting the possibility of 
Hfq assisted degradation [86]. Some sRNAs such as RhyB 
targeting sodB mRNA [87] and SgrS targeting ptsG 
mRNA primarily block translation initiation [88], though 
for many mRNAs, concomitant degradation makes 
repression irreversible [45]. In other cases such as MicC 
sRNA, degradation is the sole route of repression [89] 
(Fig. 2B). Although the vast majority of known sRNAs 
destabilize mRNA, two exceptions are the GadY and 
VR-RNA sRNAs that exert a stabilization effect on the 
gadX and colA mRNAs respectively [5, 90, 91] (Fig. 2B). 
The degradation effect exerted by natural sRNA has also 
been emulated by synthetic versions of sRNA [92]. In this 
work, Man et al. [92] developed a series of synthetic 
sRNAs designed to target 5' UTRs to trigger RNase 
E-mediated degradation and showed that the inclusion of 
natural Hfq-binding sites stabilized the synthetic sRNA in 
vivo and facilitated gene silencing. 

2.2.3 mRNA degradation control: Ribozymes 

The cleavage of mRNA catalyzed by ribozymes has also 
been shown to effect mRNA degradation. For example 
members of the glmS class combine ligand response to 
glucosamine-6-phosphate (GlcN6P) with self-cleavage 
and thus can be considered as a type of riboswitch [93] 
(Fig. 2B). The cleavage catalyzed by the glmS ribozyme 
results in rapid degradation of cleaved products [93] . Rely- 
ing on a similar mechanism, Carothers et al. [11] engi- 
neered aptazyme and ribozyme dependent cleavage of 
mRNA to control degradation for the regulation of biosyn- 
thetic pathways (Fig. 2B). 

2.2.4 mRNA degradation control: Riboswitches 

The inhibition of translation caused by riboswitches has 
been shown to indirectly induce mRNA degradation. As 
translating ribosomes can protect mRNA species from 
decay, when translation is inhibited by ligand binding to 
a riboswitch, an indirect consequence is mRNA degrada- 
tion as hypothesized for the thiM and btuB riboswitches 
[94]. However, recently the E. coli lysC riboswitch was 
shown to directly induce mRNA degradation in addition 
to inhibiting translation by ligand-dependent rearrange- 
ments resulting in exposure of an RNase E site [95]. 
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2.3 Transcriptional control of gene expression 
by RNA 

Bacterial transcriptional regulation at the RNA level 
occurs primarily through transcription termination. 
Unlike for translation and mRNA degradation, specific 
RNA sequences and structures are required to terminate 
transcription. In rfio-mediated termination, the protein 
factor iho binds to specific binding sites of mRNAs and 
traverses the mRNA in a 5'— 3' direction, eventually caus- 
ing actively transcribing polymerase to abort [96]. More 
common is intrinsic termination, where hairpin struc- 
tures placed immediately upstream of extended polyU 
sequences fold while the polymerase is paused at the 
polyU, causing polymerase to abort [97, 98]. While these 
mechanisms are ubiquitous in nature, RNA-mediated 
transcription regulation is just being recognized as a par- 
ticularly powerful point of control in synthetic biology [99] . 
Indeed, engineering RNA-mediated transcriptional regu- 
lators holds the promise of constructing RNA-only genet- 
ic networks that require no protein intermediates [99]. 

While this has only been shown as a proof-of-concept, 
potential advantages of such RNA-only networks over 
protein-based networks are faster signal propagation due 
to quick RNA degradation rates, and a smaller genetic 
footprint as regulatory RNAs can be encoded in fewer 
nucleotides than proteins. Additionally, in certain 
instances RNA-only networks possibly offer lower meta- 
bolic burden since amino acids are not consumed for 
translation; however, this has yet to be fully explored. For 
all the potential advantages, work in this area is only just 
beginning, as discussed below. 

2.3. 7 Transcriptional control: sRNA 

Transcription attenuation is a repressive sRNA-based 
mechanism where an antisense sRNA targets an attenu- 
ator sequence that lies in the 5' UTR of the target mRNA 
(Fig. 2A). The attenuator sequence can adopt two mutu- 
ally exclusive confirmations depending on the base-pair- 
ing to sRNA: in the presence of sRNA, a terminator struc- 
ture is formed in the mRNA thereby causing premature 
termination of transcription before the coding sequence 
is reached, while in the absence of sRNA, anti-terminator 
structures are formed that inhibit terminator formation 
and enable transcriptional read-through [5]. To date, this 
mechanism has been found for cis-encoded sRNAs in 
plasmids and phages of Gram-positive bacteria [5, 6] and 
has been shown to function in Gram-negative bacteria [5] . 

Of these natural mechanisms, the S. aureus plasmid 
pT181 attenuator [100] has been shown to be a versatile 
regulatory building block for RNA-based genetic networks 
[99]. In particular, independently acting, or orthogonal, 
antisense/attenuator pairs were engineered through 
mutations and shown to independently regulate GFP and 
RFP targets in the same cell. Following the example of nat- 
ural riboswitches, these orthogonal attenuators were also 



placed in tandem on the same transcript and shown to 
evaluate genetic NOR logic [99] . Importantly, these atten- 
uators were also used to construct the first example of an 
RNA-only network, a two stage transcriptional cascade, 
whereby signals were propagated directly as RNA mole- 
cules with no intermediate protein signals. In addition, the 
same aptamer fusion strategy highlighted above was used 
to produce aptamer-pT181 antisense variants that 
respond to theophylline and the MS2 coat protein, open- 
ing the door for creating sophisticated RNA circuitry that 
dynamically responds to environmental and intracellular 
cues [60] (Fig. 2A). Recently, large libraries of orthogonal 
pT181 antisense/attenuator variants were created 
through a chimeric fusion strategy to translational sRNAs, 
thus removing a critical barrier to making larger RNA- 
based genetic circuits from this system a reality [101]. 

2.3.2 Transcriptional control: Riboswitches 

Transcriptional regulation by riboswitches involves 
expression platforms that fold into terminator hairpins or 
alternative anti-terminator structures much like the tran- 
scriptional attenuators discussed above (Fig. 2A). In addi- 
tion, a class of riboswitches that regulates gene expres- 
sion via rho-dependent transcriptional termination has 
recently been discovered [102]. They are generally repres- 
sive, such as thiamin pyrophosphate (TPP) and flavin 
mononucleotide (FMN) riboswitches from Bacillus subtilis 
[103], and utilize the ligand dependent formation of ter- 
minator structures [62] (Fig. 2A). There are also a few 
known examples of activators such as the adenine and 
glycine riboswitches [66, 104] (Fig. 2A). Although ribo- 
switches are widely present in nature, the development of 
synthetic variants of riboswitches has been limited, and 
work has mainly focused on translational control due, in 
part to the limited structural characterization of tran- 
scriptional expression platforms. Recently, however, 
Wachsmuth et al. [13] used an in silico guided approach 
to screen variants of riboswitch designs that varied in a 
spacer sequence between the aptamer and a synthetic 
transcription expression platform. Several variants in the 
spacer sequence were shown to functionally activate 
transcription in E. coli in the presence of the theophylline 
ligand. Another approach harnessed the modular nature 
of aptamers and expression platforms. Ceres et al. [105] 
showed that naturally occurring and synthetic aptamers 
could be effectively fused with three different naturally 
occurring expression modules of riboswitches from 
B. subtilis to create chimeric transcriptional riboswitches 
that respond to tetracycline and theophylline. 

2. 3. 3 Transcriptional control: Leader peptides 

Leader peptides are naturally occurring cis-regulatory ele- 
ments that couple translation of a small coding sequence 
to the transcription of downstream genes (Fig. 2A). Two 
notable examples of natural leader peptide regulators 
integrate tryptophan-dependent translation of a leader 
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peptide to the control of downstream genes in the tip and 
tna operons. In the tip operon, a 14 amino acid leader pep- 
tide coding sequence contains two tandem tryptophan 
codons, and is directly upstream of a sequence that can 
form either a transcriptional terminator or an anti-termi- 
nator based on the concentration of charged tryptophan 
tRNAs [106] (Fig. 2A). In contrast, the tna operon relies 
upon ribosome stalling in a leader peptide to prevent 
transcriptional termination by physical occlusion of a Rho 
termination factor-binding site (mt site) [106] (Fig. 2A). 
Recently, these mechanisms have been engineered to 
sense the presence of unnatural amino acids, and to 
expand the repertoire of sRNA-mediated transcriptional 
regulators available to synthetic biologists [107, 108]. In 
the first example, Lui et al. [107] designed variants of the 
leader peptide sequence to contain "blank codons" that 
do not naturally encode amino acids, thereby inducing 
ribosomal stalling in both the tip and tna leader peptides 
(Fig. 2A). By introducing tRNAs to decode these blank 
codons, and aminoacyl-tRNA-synthetases to charge 
these tRNAs with unnatural amino acids, they showed 
that they could couple the presence of unnatural amino 
acids to either transcriptional activation (tna leader pep- 
tide) or repression (tip leader peptide). In the second 
example, Liu et al. [108] demonstrated that sRNAs 
could be used to regulate the translation of the leader pep- 
tides themselves, thus converting this translational regu- 
lation into transcriptional regulation of the downstream 
sequence (Fig. 2A). 

2.3.4 Transcriptional control: CRISPRs systems 

While natural CRISPR systems are involved in complex 
bacterial defense systems, recently they have been engi- 
neered to provide a powerful mechanism of RNA-meditat- 
ed transcription control. In an exciting new advance, Qi et 
al. [109] describe what they call CRISPR interference 
(CRISPRi), a system engineered from the type II CRISPR of 
Stieptococcus pyogenes (Fig. 2A). In this system, the Cas9 
protein forms a complex with a small guide RNA (sgRNA) 
that contains a region complementary to a target genom- 
ic DNA sequence. In the native system, the Cas9:sgRNA 
complex binds to the target DNA and the nuclease activi- 
ty of Cas9 then cleaves the DNA. However, in this case, the 
authors used a catalytically dead Cas9 mutant (dCas9), 
which still binds to the target and was shown to block RNA 
polymerase (RNAP) binding during either transcription 
initiation or elongation. This work effectively opens anoth- 
er route for creating RNA-only genetic circuitry by config- 
uring networks of sgRNAs that control their own tran- 
scription in sophisticated patterns, and could be synergis- 
tic with the sRNA-based networks described above. 

2.4 Regulation of protein function by RNA 

In addition to the regulation of gene expression, RNA is 
able to regulate post-translationally by directly interact- 



ing with proteins. Modulation of protein function is typi- 
cally achieved by sequestration: for example, the 
E. coli a 70 RNAP holoenzyme can be sequestered by the 
6S RNA, which is thought to mimic the conformation of 
a 70 DNA promoters elements [7]. Similarly the CsrA pro- 
tein involved in mRNA translation and degradation can 
be sequestered by the CsrB RNA [7]. 

RNA-protein interactions have successfully been 
engineered using the ability of RNA aptamers to recog- 
nize and bind a variety of proteins. For example, aptamers 
able to bind the TetR transcriptional regulator have been 
successfully generated and placed in the 5' UTR of 
mRNAs to confer post-transcriptional regulation of gene 
expression by TetR [110]. Alternatively, these aptamer- 
protein interactions have been harnessed to build one and 
two-dimensional RNA scaffolds in vivo that include bind- 
ing sites for tagged proteins, thus enabling intracellular 
localization of multiple proteins [111]. These RNA scaf- 
folds were used to organize and enhance the production 
of a hydrogen-producing enzymatic pathway [111]. 

3 Tools for engineering the RNA sequence- 
structure-function relationship 

It is clear that the rich versatility of RNA gene regulatory 
function is intimately tied to RNA structure. In fact, it is 
precisely the ability of single stranded RNA molecules to 
fold back on themselves with strong base-pairing interac- 
tions that allows fundamental processes of gene expres- 
sion to be blocked or allowed as discussed above. In addi- 
tion, unpaired nucleotides can fold into complex three- 
dimensional topologies that create specific binding sites 
for other cellular RNAs, small molecules, and proteins 
that, when bound, can alter the underlying fold and regu- 
latory function of the RNA [4] . Thus our capacity to engi- 
neer with RNA boils down to our ability to engineer its 
fundamental sequence-structure-function relationship. 

Fortunately, we have a rapidly emerging understand- 
ing of this relationship based on a strong biophysical 
foundation that has already led to the beginnings of a sys- 
tem of computational RNA design. Much like for proteins, 
a given RNA structure is divided into secondary and ter- 
tiary structures, the latter defined as the three-dimen- 
sional arrangement of its atoms [112]. Unlike for proteins, 
however, secondary structures often form faster and with 
more stable free energies than tertiary structures [113]. 
This critical observation has allowed the decoupling of 
the RNA secondary and tertiary structure folding prob- 
lems, and great progress has been made in predicting the 
ensemble of secondary structures of an RNA sequence 
using experimentally determined thermodynamic base- 
pairing parameters [112, 114]. 

Since RNA-based gene regulation often involves the 
conditional formation of specific base-pairing patterns, 
RNA secondary structure prediction algorithms can be 
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immediately converted into computational design tools 
for engineering RNA-based gene regulators. In addition, 
recent advances in biochemically characterizing RNA 
structures and interactions have created experimental 
approaches that can improve structural modeling and 
enable observation of ligand-mediated structural switch- 
ing, but these advances are only just beginning to be used 
for RNA engineering. In this section we give an overview 
of these computational and experimental methods in the 
context of engineering RNA-based gene regulation. 

3.1 Computational methods for predicting 
RNA structure 

The most well established methods for computationally 
predicting RNA structures focus on calculating the sec- 
ondary structure partition function for an RNA sequence. 
The partition function, defined from statistical thermody- 
namics as the Boltzmann-weighted sum over all possible 
RNA folds, mathematically characterizes the ensemble of 
structures in which a given RNA sequence can exist in 
solution. Once calculated, the partition function is the key 
to equilibrium structural properties of the RNA molecule 
and can be used to determine the minimum free energy 
(MFE) RNA structure, suboptimal folds that are close to 
the MFE, the probability that a given region is base paired 
over the ensemble of folds [112], and can even be used to 
statistically sample folds that occur in the ensemble [115]. 
There are several algorithms that perform these calcula- 
tions that differ mainly in implementation details and 
thermodynamic parameter sets. The most commonly 
used algorithms that also have convenient web servers 
include RNAStructure [116], ViennaRNA [117], and 
Unafold [114]. 

The ability to model RNA secondary structures has 
been applied to allow more rational engineering of RNA 
regulators. Most recently, Rodrigo et al. [12] used a search 
algorithm based on ViennaRNA to assess potential 
designs of taRNA/crRNA riboregulator pairs as discussed 
previously. By searching through sequences and model- 
ing their folding and interactions, they were able to gen- 
erate six orthogonal taRNA/crRNA pairs that showed up 
to 11 -fold activation of translation. Similarly, in order to 
reduce the number of constructs that had to be experi- 
mentally tested, Qi et al. [42] used the M-fold algorithm 
[118] to assess the folding of designed aptamer and non- 
coding RNA (ncRNA) fusions to find versions that were 
predicted to interfere with the structure of the ncRNA in 
the absence of ligand [60]. Additionally, rational design 
algorithms have used folding algorithms to develop syn- 
thetic riboswitches that activate transcription [13]. In this 
work, Wachsmuth et al. [13] computationally generated 
designs that varied within a spacer and a expression plat- 
form sequence, while maintaining aptamer sequence, 
and assessed the MFE structure based on several criteria 
designed to measure functionality. Finally, folding algo- 



rithms have been used to develop design tools focused on 
the role of general RNA structure in gene regulation. In 
particular, Salis et al. [10] developed the RBS calculator, an 
algorithm that uses the Nupack suite of RNA folding algo- 
rithms [119] to calculate free energies of competing RNA 
structures that could prevent ribosome binding to the 
RBS. Using this information, the RBS calculator models 
the effective translation initiation rate of a given mRNA 
sequence competing for ribosome binding within a pool 
of cellular mRNAs [10]. The RBS calculator was shown to 
predict measured strength of GFP expression from a 
library of 5' UTRs within a factor of 2.3 over a 100 000-fold 
range and could be used in a forward design mode to 
design an RBS within an mRNA context to produce a 
specified amount of protein expression [10]. 

While excellent progress has been made in using com- 
putational RNA structure prediction tools to engineer 
RNA gene regulatory mechanisms, we emphasize that 
their inherent approximations should be kept in mind 
when applying them. First and foremost, most of these 
algorithms only take into account secondary structures 
formed from canonical base-pairing interactions, 
although programs such as MC-Fold and MC-Sim have 
begun to address this problem [120]. Even then, most 
algorithms do not predict the formation of pseudoknots, 
which occur when nucleotides in a loop pair with a region 
outside the local helices that close the loop [121]. 
Although a new wave of algorithms such as Turboknot 
[122] and HotKnots [123] allow pseudoknot formation, 
these algorithms still do not incorporate the rich set of 
non-canonical base-pairing interactions known to form 
functional motifs in a wide array of natural RNAs [124]. 
The prediction and design of these motifs remain at the 
forefront of computational RNA structure methods, so 
such motifs are not currently designable using these tools. 

While these approximations should be kept in mind 
when using secondary structure folding algorithms to 
engineer RNA folds, it is also important to realize that they 
model folding conditions where the RNAs have reached 
thermodynamic equilibrium. Other frontiers in RNA 
structure-prediction include taking into account co-tran- 
scriptional folding kinetics. Co-transcriptional folding has 
been shown to dynamically induce structural rearrange- 
ments in the mRNA being transcribed and is particularly 
important for structural switches such as riboswitches, 
ribozymes, and leader peptide attenuators [125]. Compu- 
tational methods such as Kinefold have begun to provide 
tools to study RNA folding pathways during transcription, 
albeit with some of the same inherent limitations regard- 
ing approximate thermodynamic parameters and non- 
canonical base pairs discussed above [126]. Recently, 
Kinefold was used in an RNA engineering context to 
design appropriate sequence contexts for proper folding 
of ribozymes and aptazymes engineered to control the 
mRNA degradation rate, and thus enzyme levels, in a 
biosynthetic pathway [11]. We anticipate increased appli- 
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cation of co-transcriptional folding algorithms to design- 
ing RNA folds, particularly those that dynamically change 
state to control gene expression. 

3.2 Experimental methods for characterization 
of mRNA structure 

Despite the many advances in computational structure 
methods, the approximations discussed above mean that 
predicted structures are only a model of what is actually 
present inside the cell. Although numerous RNA species 
such as riboswitches [127, 128] and ribozymes have been 
resolved structurally by X-ray crystallography and NMR, 
due to the low-throughput nature and the specialization 
required for these methods, the number of structures is 
limited, and to date there are no crystal structures of 
sRNAs. Experimental chemical and enzymatic probing 
methods offer an alternate way to infer RNA structures 
and interactions with other molecules, and the potential 
to couple these methods with next-generation sequenc- 
ing (NGS) opens the door to ultra high-throughput RNA 
structural determination (see Section 4). 

Experimental RNA structure probing relies on the use 
of enzymes or chemicals that either cleave or chemically 
modify RNA in a structure-dependent fashion (Fig. 3). 
A wide variety of these chemical and enzymatic probes 
exist, each with different selectivities for a particular RNA 
sequence or structure [129, 130]. These techniques are 
commonly used on RNAs that have been folded in vitro, 
the structure of which is determined by quantifying the 
number of modification or cleavage events at each 
nucleotide, which is inherently linked to the underlying 
structure of the RNA. Modified or cleaved positions can 



be identified directly by gel electrophoresis or through 
reverse transcription, which is blocked by modification or 
cleavage [131]. Reverse transcription creates a pool of 
truncated DNA species whose 3' ends marks the location 
of a chemical modification or cleavage [132]. These 
reverse transcription products can be read off using a 
variety of techniques, and there has been considerable 
effort in making chemical probing experimental protocols 
relatively straightforward and high-throughput by using 
capillary electrophoresis [133] and even NGS [17, 18, 20]. 

Once reverse transcription products have been quan- 
tified, there are now robust methods of converting the dis- 
tributions of modification position into an estimate of the 
"reactivity" of each nucleotide to the chemical probe [19, 
134]. Reactivities can then be used to infer RNA struc- 
tures and RNA-ligand interactions. Recent progress in 
this area has utilized experimental reactivities as con- 
straints in RNA-folding algorithms by first converting 
them into pseudo-free energies that are incorporated 
along with the thermodynamic free energy parameters 
[135]. The implementation of this approach in the 
RNAStructure package has been shown to increase the 
accuracy of secondary structure predictions [135], and 
work is being done to establish best practices when uti- 
lizing and interpreting this information [136, 137]. Addi- 
tionally, reactivity constraints have now been included in 
versions of the RNAStructure 's pseudoknot prediction 
tool [121] as well as the SeqFold structure predication 
algorithm [138]. 

An alternative to chemical and enzymatic probing is 
in-line probing, a technique that harnesses the ability of 
RNA to non-enzymatically cleave itself at differential rates 
according to the RNA structure [139]. RNA molecules 



1) Chemically Modify 



4) Infer Structure 




2) Convert to DNA 
and Sequence 



3) Analyze 



Figure 3. Chemical probing of RNA 
structure. Specific chemicals covalently 
modify folded RNAs at unpaired and 
flexible nucleotide positions. These mod- 
ified RNAs are converted into DNA via 
reverse transcription, which is blocked 
by the modifications, thus creating trun- 
cated single-stranded DNA (ssDNA) 
species. These ssDNA species are then 
sequenced, and locations of chemical 
modifications identified to derive a 
measure of nucleotide "reactivity" to the 
chemical. Reactivities are then used to 
infer structural properties of the input 
RNAs. 
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such as riboswitches are transcribed and folded in vitro 
and incubated in the presence or absence of ligand for 
around 48 h. Cleaved products can then be analyzed 
directly by gel electrophoresis and the structure elicited 
[139]. 

Although the application of experimental structure 
probing to RNA engineering is just beginning, the current 
research does highlight the enormous potential for struc- 
tural and functional insight. First of all, these experiments 
provide validation and improvements of computational 
predictions, thus increasing the accuracy of structural 
information forming the basis for rational design. For 
example, Liu et al. [108], in their design of tna leader pep- 
tide adaptors, used experimental structure determination 
to provide structural validation that secondary structures 
of the wild-type system could be readily disrupted and 
still produce a functional regulator. Secondly, experimen- 
tal structure probing is sensitive to more than just sec- 
ondary structure, and can also be used to analyze RNA- 
ligand interactions. In fact, one common chemical prob- 
ing technique called selective 2'-hydroxyl acylation ana- 
lyzed by primer extension (SHAPE) produces lower 
reactivities at positions that participate in RNA-ligand 
interactions, and has been used to identify small-mole- 
cule-aptamer interactions in riboswitches [140] and RNA- 
protein interactions [141]. In a recent use of this feature in 
an engineering context, Qi et al. [60] showed how SHAPE 
probing could give structural mechanistic validation of 
the engineered IS10 RNA-OUT-theophylline aptamer 
fusion. 

Perhaps the two most promising aspects of structural 
probing for RNA engineering are the adaption of these 
techniques to high-throughput sequencing platforms and 
the improvements in the chemistry and range of applica- 
tion of the probes. In particular, recent innovations in NGS 
have enabled the multiplexing of chemical probing exper- 
iments, thus significantly increasing the throughput of 
structural probing experiments. These high-throughput 
approaches include parallel analysis of RNA structure 
(PARS) [17] fragmentation sequencing (FragSeq) [18], and 
selective 2'-hydroxyl acylation analyzed by primer exten- 
sion sequencing (SHAPE-Seq) [20] and are discussed in 
detail in Section 4. In addition, newer generations of 
SHAPE probes such as 2-methylnicotinic acid imida- 
zolide (NAI) and 2-methyl-3-furoic acid imidazolide (FAI) 
allow reactions to be carried out at physiological condi- 
tions and in vivo [142]. This opens the door to character- 
ize engineered RNA regulators within the myriad of 
potential interactions in the cellular milieu. 

Thus chemical probing has the power to not only 
improve the accuracy of structural modeling, but also to 
aid in our understanding of more complex RNA interac- 
tions not accessible through computation alone. Experi- 
mental structural probing can help elucidate key mecha- 
nisms of ligand binding, intermolecular interaction, and 
conformational switching that lie at the heart of engi- 



neering complex modes of RNA regulation. We anticipate 
these technologies to be instrumental in future RNA engi- 
neering efforts. 

4 Characterizing RNA gene regulators 
with next-generation sequencing 

NGS refers to a suite of technologies that have created an 
explosion in the throughput of DNA sequencing (for a 
review of the technology see [25]). At their core is a 
sequencing-by-synthesis process that can be spatially 
arrayed and imaged for millions of molecules at a time. In 
fact, the high-throughput of NGS has led to the emer- 
gence of a whole family of creative applications, called 
"sequence census" methods, that use the capability to 
count DNA molecules to ask questions about biological 
systems [143]. 

The principle of sequence census methods is simple: 
(i) use biochemical processes to convert a question of 
interest into a pool of DNA molecules, (ii) "count" this pool 
with NGS, (iii) statistically correct for any biases in the 
conversion and sequencing, and (iv) analyze the cor- 
rected counts to answer the question. For example, if the 
question is "what RNAs are being expressed in a cell pop- 
ulation?" then the biochemical conversion would consist 
of extracting RNA from the population and using reverse 
transcription to convert it into DNA, while statistical 
analysis might include normalizing counts to internal ref- 
erences. In fact, such a technique, called RNA-Seq, is well 
established [15, 16] and similar techniques are applicable 
to measuring RNA structures and interactions with a 
wide array of important cellular molecules (Fig. 4). 

4.1 Characterizing RNA abundance 

RNA-Seq uses NGS to measure RNA abundance [15, 16] 
and involves isolation of an RNA population, fragmenta- 
tion of RNA, and conversion into DNA by random prim- 
ing. After platform-specific library preparation steps, the 
library of fragments is then sequenced using the chosen 
NGS platform. The experimental steps are well estab- 
lished in a number of vendor-specific and published pro- 
tocols [144], and there are robust computational pipelines 
available for analyzing read counts to quantify RNA abun- 
dance and differential gene expression among different 
conditions [145]. 

RNA-Seq is an obvious and powerful component of 
the synthetic biology characterization toolbox, especially 
for engineered RNA gene regulators. Since all cellular 
RNAs are examined simultaneously, RNA-Seq has the 
potential to characterize both introduced synthetic 
genetic circuitry and host gene expression patterns, far 
beyond what the typically used fluorescent protein 
reporters allow. This capability will allow synthetic biolo- 
gists to monitor and even engineer interactions with host 
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Figure 4. Next-generation sequencing can characterize RNA abundances, structures and interactions across the cell. Names of established NGS-tech- 
niques that characterize these features are highlighted next to the interactions examined. 



networks, which could become critical as the field moves 
to engineer larger and more sophisticated genetic net- 
works. In addition, RNA-Seq has been used to character- 
ize RNA abundance pools generated by SELEX to enable 
more efficient generation of RNA aptamers [146]. 

4.2 Characterizing RNA structures 

Beyond measuring the abundance of a particular RNA 
sequence, NGS is providing a platform to generate struc- 
tural information. As mentioned above, there currently 
exist three distinct approaches to the high-throughput 
characterization of RNA structure; PARS [17, 147] 
FragSeq [18], and SHAPE-Seq [19, 20, 148]. PARS and 
FragSeq are similar in that they use structure dependent 
nuclease cleavage to create a distribution of RNA frag- 
ments that are sequenced with a variant of the RNA-Seq 
protocol. By comparing fragment distributions generated 
from nucleases with different structural sensitivity, struc- 
tural signatures can be discerned for thousands of RNA 



species in vitro. In contrast, SHAPE-Seq utilizes a more 
sensitive chemical probe, influenced by both secondary 
and tertiary RNA structure, to first chemically modify 
RNAs in vitro, followed by conversion to cDNA and 
sequencing as described above [19, 20, 148]. Further- 
more, the SHAPE-Seq protocol and data analysis method 
[148] should be generally applicable to other chemical 
probes that interrogate other aspects of RNA structure 
such as solvent accessibility [149] and in vivo structures. 
Despite these differences, it should be mentioned that 
each of these techniques are in their infancy, and we ful- 
ly expect dramatic improvements in accuracy, applicabil- 
ity and throughput. 

4.3 Characterizing RNA transcription 
and translation 

Ribosomal profiling is an NGS-based experiment that 
determines the number of ribosomes translating a partic- 
ular mRNA species in vivo [21 ] . The technique is based on 
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capturing actively translating ribosomes with affinity 
pull-downs, using exonucleases to chew up RNA 
sequences not bound by protein, then releasing ribo- 
some-bound sequences for subsequent library prepara- 
tion and sequencing (for a protocol see [150]). Analyzing 
where ribosome-bound sequences fall on each mRNA in 
the transcriptome allows the quantification of the relative 
translation rate on an mRNA species at a given time [151], 
and also allows researchers to characterize sites of non- 
canonical translation [152]. Similarly, the position and ori- 
entation of RNAP and therefore the effective transcription 
rate can be determined by global run-on sequencing 
(GRO-Seq) [22], an analogous technique to ribosomal pro- 
filing that targets actively transcribing RNAPs. In an engi- 
neering context, profiling ribosome and RNAP positions 
could be a powerful tool for synthetic biology. In particu- 
lar, these techniques can be used to understand and trou- 
bleshoot network dynamics by quantifying the translation 
and transcription rate of individual nodes in a network. 
Moreover, by comparison to endogenous RNA species 
they can help assess the burden of an implemented net- 
work on the population of RNAP and ribosomes, giving 
complementary information to RNA-Seq experiments. 
Finally, these methods can be used to functionally char- 
acterize RNA regulators of gene expression by measuring 
the effect of RNA regulators on transcription and transla- 
tion rates of target mRNAs. 

4.4 Characterizing RNA interactions 

Much like for ribosomal profiling, RNA's interaction with 
other cellular proteins can be determined using NGS. 
RNA immunoprecipitation (RIP-Seq) uses immunoprecip- 
itation of an RNA-protein complex and subsequent 
reverse transcription and sequencing to identify which 
RNA sequences interact with a protein of interest [23]. In 
addition, intermolecular and intramolecular RNA-RNA 
interactions can be resolved using cross-linking, ligation, 
and sequencing of hybrids (CLASH) [24]. The CLASH 
method relies upon purification of RNA-RNA hybrids 
mediated by cross-linking to a protein-binding partner. 
These RNA-protein complexes are then immunoprecipi- 
tated and the RNA-RNA hybrids ligated, reverse tran- 
scribed and sequenced. Because this technique relies 
upon a protein interaction to the RNA-RNA hybrid it has 
so far only been used to study ribonucleoprotein com- 
plexes. Expansion of this technique to enable generic 
purification of RNA-RNA hybrids could dramatically 
increase its potential use for synthetic biology. 

In summary, NGS provides an array of high-through- 
put tools to study RNA within the context of the cell. 
Although its application to RNA engineering is yet to 
be fully realized, we believe it is an immensely diverse 
tool that has the power to fundamentally change the 
approach to characterizing RNA-based gene regulatory 
systems. 
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5 Conclusion 

We reviewed and presented the case for RNA's powerful 
role in the engineering of biological systems. We believe 
that the combination of being a versatile, designable, and 
now broadly characterizable regulator of gene expression 
places RNA in a unique and central position as a master 
substrate for engineering gene expression. Although 
much work remains to be done to link these three power- 
ful features and realize RNA's potential, as reviewed 
above, important strides have been made. While conceiv- 
ably RNA synthetic biology is still mainly in the phase of 
engineering basic regulatory parts, its emerging transi- 
tion toward engineering networks and higher-order regu- 
latory systems holds great promise for synthetic biology 
as a whole. 
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Biotechnology Journal's latest Special Issue on "Biomolecular Engineering" is based 
peer-reviewed papers derived from some of the best presentations at the 4 th Inter- 
national Conference on Biomolecular Engineering. The cover is a snapshot of solvated 
HER4/ErbB4 kinase domain with juxtamembrane region; image provided by Ravi 
Radhakrishnan. 
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